X-Ray Scaling Relations of Galaxy Groups in a Hydrodynamic 

Cosmological Simulation 



Romeel Dave ^ 
Steward Observatory, University of Arizona, Tucson, AZ 85721 

Neal Katz 

Astronomy Department, University of Massachusetts, Amherst, MA 01003 

and 

David H. Weinberg 

Astronomy Department, Ohio State University, Columbus, OH 43210 



ABSTRACT 

We examine the scalings of X-ray luminosity, temperature, and dark matter or galaxy velocity 
dispersion for galaxy groups in a ACDM cosmological simulation, which incorporates gravity, gas 
dynamics, radiative cooling, and star formation, but no substantial non-gravitational heating. In 
agreement with observations, the simulated Lx — c and Lx — Tx relations are steeper than those 
predicted by adiabatic simulations or self-similar models, with Lx oc cr*'^ and Lx oc T^^ for mas- 
sive groups and significantly steeper relations below a break at cr « 180 km/s {Tx « 0.7 keV). 
The Tx — u relation is fairly close to the self-similar scaling relation, with Tx oc ct^'^^, provided 
that the velocity dispersion is estimated from the dark matter or from > 10 galaxies. The en- 
tropy of hot gas in low mass groups is higher than predicted by self-similar scaling or adiabatic 
simulations, and it agrees with observational data that suggest an "entropy floor." The steeper 
scalings of the luminosity relations are driven by radiative cooling, which reduces the hot (X-ray 
emitting) gas fraction from 50% of the total baryons at cr « 500 km/s to 20% at a ~ 100 km/s. 
A secondary effect is that hot gas in smaller systems is less clumpy, further driving down Lx- 
A smaller volume simulation with eight times higher mass resolution predicts nearly identical 
X-ray luminosities at a given group mass, demonstrating the insensitivity of the predicted scal- 
ing relations to numerical resolution. The higher resolution simulation predicts higher hot gas 
fractions at a given group mass, and these predicted fractions are in excellent agreement with 
available observations. There remain some quantitative discrepancies: the predicted mass scale 
of the Lx — Tx and Lx — cr breaks is somewhat too low, and the luminosity-weighted tempera- 
tures are too high at a given a, probably because our simulated temperature profiles are flat or 
rising towards small radii while observed profiles decline at r < 0.2i?vir- We conclude that radia- 
tive cooling has an important quantitative impact on group X-ray properties and can account for 
many of the observed trends that have been interpreted as evidence for non-gravitational heating. 
Improved simulations and observations are needed to understand the remaining discrepancies and 
to decide the relative importance of cooling and non-gravitational heating in determining X-ray 
scalings. 

Subject headings: cosmology: observations — cosmology: theory — X-rays: galaxies: clusters — galaxies: 
clusters: general 
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1. Introduction 

Galaxy groups, i.e. bound systems having typ- 
ically a few Li, galaxies, contain the majority of 
galaxies in the Universe, and hence seem crucial 
to understanding the processes of galaxy forma- 
tion. However, observational and theoretical stud- 
ies have tended to focus either on the smaller scale 
of individual galaxies or on the larger scale of rich 
clusters. Groups present an observational chal- 
lenge, since projection effects make them difficult 
to identify unambiguously in optical imaging sur- 
veys and their X-ray emission is faint, emerging 
predominantly at lower energies that are contam- 
inated by foreground Galactic emission. Sensi- 
tive X-ray instruments and group catalogs created 
from large redshift surveys are transforming this 
observational situation, and they have revealed 
puzzling departures in the group mass regime from 
the scaling relations obeyed by rich clusters. In 
this paper we use a large hydrodynamic cosmo- 
logical simulation to examine the predicted X-ray 
scaling relations of galaxy groups, reaching to the 
mass regime of poor clusters. We adopt an infla- 
tionary cold dark matter scenario with a cosmo- 
logical constant (ACDM). 

The canonical approach to understanding 
groups treats them as scaled down versions of rich 
clusters. The simplest model of a cluster, in turn, 
is a sphere of Virial-temperature gas punctuated 
by old galaxies, with perhaps a cooling flow onto 
the central cD galaxy. In this model, the gas cool- 
ing time is longer than a Hubble time everywhere 
except near the center, and thus it is predicted 
that intracluster and intragroup media should fol- 
low "self-similar" scaling relations: Tx oc cr|g^j, 
Lx oc T"^. and Lx oc cTgai, where Lx is the to- 
tal X-ray luminosity, Tx is the gas temperature 
(presumed to be the halo Virial temperature), and 
(Tgai is the velocity dispersion of cluster galaxies 
(presumed to trace the system's total mass). The 
first relation arises from equating the gas thermal 
energy to the galaxies' kinetic energy (i.e. hy- 
drostatic equilibrium), the second from combining 

1/2 

Lx oc MT^ for free-free emission (where M is 
the system mass), M oc (t|,j from the Virial the- 
orem, and the Tx — Cged relation, and the third 
from combining the previous two scaling relations 

^Hubble Fellow, rad@as.arizona.edu 



(sec e.g. Navarro, Frenk, & White 1995). 

While observed properties of the most massive 
clusters follow these relations (Allen & Fabian 
1998; Xu, Jin & Wu 2001), the self-similar model 
breaks down as one progresses to smaller sys- 
tems, with observations indicating lower than 
expected luminosities for a given temperature 
or velocity dispersion. For 207 clusters ob- 
served with Einstein, White, Jones, & Forman 
(1997) found Lx oc T^. Extending such obser- 
vations to groups has yielded conflicting results; 
Mulchaey & Zabludoff (1998, hereafter MZ98) 
found Lx oc T^^, consistent with clusters, while 
Helsdon & Ponman (2000, hereafter HPOO) and 
Xue & Wu (2000) found much steeper relations, 

r ^ m4.9-5.6 

Lx oc I-^ 

In an influential paper, Ponman, Cannon, 
& Navarro (1999, hereafter PCN99) calcu- 
lated the specific entropy of gas at 10% of the 
Virial radius, and found an "entropy floor" of 
~ lQOh~^^^ keV cm~^ among their smallest sys- 
tems. This entropy floor reduces the X-ray lumi- 
nosity by lowering the gas density progressively 
more in smaller systems, and hence reproduces 
the steep Lx — Tx relations seen in groups. They 
suggested that the entropy floor arises from en- 
ergy injection by supernova winds, which might 
concurrently explain the enrichment level of the in- 
tracluster medium (ICM). Though there arc large 
uncertainties for the handful of groups that de- 
marcate the entropy floor, this result nevertheless 
spawned numerous papers investigating the possi- 
ble physical sources of "pre-heating," constraints 
on parameters of the pre-heating model such as en- 
ergy per baryon and energy injection redshift, and 
implications of pre-heating for other group and 
cluster observablcs. For instance, Lloyd-Davies, 
Ponman, & Cannon (2000) concluded that a tem- 
perature increase of 0.3 keV per baryon is neces- 
sary to explain the observed entropy floor, while 
Tozzi & Norman (2001) considered a similar uni- 
form energy injection model and determined that 
the heat injection could be w 0.1 — 1 keV, depend- 
ing on when and how it was injected. 

Uniform heat injection across all baryons is 
ruled out by temperature measurements of the 
Lya forest; these indicate T ss 2xl0~^ keV 
for most of the baryons at z ~ 3 (Schaye et al. 
1999; Bryan 2000), and a significant fraction of 
baryons at ^ are even cooler (Ricotti, Gnedin, 
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& ShuU 2000; Dave & Tripp 2001). Thus, 
heat injection must be confined in and around ha- 
los. Semi-analytic models of this scenario require 
much greater energy injection, typically 1 — 3 kcV 
per baryon within halos (Wu, Fabian, & Nulsen 
2000; Bower et al. 2001) or entropy injection of 
~ 300 - 400 kcV cm-2 (Babul et al. 2002), to 
reproduce the PCN99 entropy fioor, a result sup- 
ported by adiabatic simulations of clusters with 
pre- heating (Bialck, Evrard, & Mohr 2001; Bor- 
ganietal. 2001). Supernovae seem unable to pro- 
vide this much energy, so another energy source 
such as AGN-driven galactic winds is necessary 
(Valageas & Silk 1999; Kravtsov & Yepes 2000; 
Wu, Fabian, & Nulsen 2000; Bower et al. 2001; 
Pipino et al. 2002) . While there is preliminary ev- 
idence for strong winds from Lyman break galax- 
ies (Pettini et al. 2001), believed to evolve into 
cluster galaxies (Governato et al. 2001), the en- 
ergy they inject seems unlikely to be as high as 

1 keV per baryon. The pre-heating model has 
also been invoked to explain the faintness of the 
soft X-ray background (Pen 1999; Wu, Fabian, & 
Nulsen 2001), though numerical simulations in- 
dicate that these constraints can probably be sat- 
isfied without non-gravitational heating (Dave et 
al. 2001; Phillips, Ostriker, k Cen 2001; Croft et 
al. 2001). In short, the pre-heating model is able 
to reproduce group and cluster scaling relations, 
but only at the expense of invoking a mysterious 
energy source that docs not naturally arise from 
any firmly established physical process of galaxy 
formation (though see Babul et al. 2002, for some 
interesting hypotheses) . This is the essence of the 
"ICM energy crisis" described by Tozzi (2001). 

A fundamental assumption underlying the pre- 
heating model is that groups, in the absence of 
such heating, would be self-similarly scaled-down 
versions of large clusters. If some process inher- 
ent in galaxy formation breaks self-similarity, it 
may be that pre-heating is not required. Grow- 
ing evidence suggests that radiative cooling could 
be this process. Bryan (2000) used an ana- 
lytic model together with the observed depen- 
dence of stellar mass fraction on cluster temper- 
ature to argue that the observed scaling relations 
can be explained without non-gravitational heat- 
ing. Voit & Bryan (2001) further argued that 
cooling removes low-entropy gas during the for- 
mation of cluster galaxies, leaving an effective en- 



tropy floor of ^ 100 koV cm^^ for groups (Voit 
et al. 2002). Hydrodynamic simulations sup- 
port these analytic arguments. Muanwong et al. 
(2001) found that simulations without cooling fol- 
lowed the self-similar relations but that the addi- 
tion of cooling steepened and lowered the Lx—Tx 
relation, bringing it into much better agreement 
with observations. Thomas et al. (2002) showed 
that adiabatic simulations underpredict the ampli- 
tude of the Tx — <y relation observed recently with 
Chandra but that the addition of either radiative 
cooling or pre-heating could resolve the discrep- 
ancy. Lewis ct al. (2000) simulated a Virgo-sized 
cluster with and without cooling and found that 
the inclusion of cooling and subsequent star forma- 
tion had an important impact on the X-ray emit- 
ting gas throughout the entire cluster, not just in 
the central region where the present-day cooling 
time is short. These results clearly point towards 
cooling as a non-negligible process in group and 
cluster formation, and they motivate a more thor- 
ough investigation of cooling as an alternative to 
pre-heating for explaining the observed scaling re- 
lations. 

In this paper we analyze a cosmological hy- 
drodynamic simulation that incorporates cooling, 
star formation, and (weak) feedback. Our sim- 
ulation has enough dynamic range to probe the 
mass range from poor groups to small clusters 
(cr « 100 — 550 km/s), while resolving sub-i* 
galaxies within a random cosmological volume. 
We describe the simulation and analysis proce- 
dures in §2. We present our simulated scaling 
relations in §3 and discuss the physical origin of 
the departures from self-similarity. We examine 
the average profiles of various physical quantities 
in our groups in §4. We go on to consider the 
possible effects of metallicity, velocity bias, and 
surface brightness thresholds in §5. These effects 
introduce some uncertainty in our predictions of 
observable scaling relations; we discuss our "best 
guess" predictions and compare them to obser- 
vations in §5.4. Section 6 presents our conclu- 
sions. In summary, we find that radiative cooling 
— in particular, the greater efficiency of cooling 
in smaller systems — can produce departures from 
self-similar scaling that quantitatively mimic those 
of a pre-heating model. Thus, observed scaling re- 
lations do not necessarily imply that groups have 
experienced substantial non-gravitational heating 
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or entropy injection. 

2. Simulation and Group Identification 

We simulate a random volume with 

a 7/i~"^kpc (equivalent Plummer, comoving) grav- 
itational softening length, assuming a ACDM cos- 
mological model with fl„i = 0.4, fl\ = 0.6, 
h = ffo/100 kms-i Mpc"^ = 0.65, n = 0.95, 
and erg = 0.8. We employ 144^ dark matter 
and 144^ gas particles, yielding particle masses of 

TObary = 8.5 X IO^Mq and mdark = 7.2 X IO^Mq. 

The input physics includes cooling, star forma- 
tion, and thermal feedback. Because the feedback 
energy is deposited locally in the dense gas giv- 
ing rise to star formation, it is usually radiated 
away quickly and produces little long-term heat- 
ing of the surrounding medium. Other than this 
feedback, we do not inject any non-gravitational 
heat or entropy. The simulation was evolved from 
^ = 49 ^ using Parallel TrecSPH (Dave et al. 
1997). A detailed discussion of the algorithms, in- 
cluding the star formation and feedback prescrip- 
tions, is given by Katz, Weinberg, & Hcrnquist 
(1996). Relative to the simulation of Muanwong 
et al. (2001), who investigate similar issues with 
a similar numerical approach, our simulation has 
a mass resolution that is higher by about a fac- 
tor of four and gravitational force resolution that 
is higher by a factor of 4-7 (depending on red- 
shift), but a simulation volume smaller by a factor 
of eight. Our simulation is therefore better suited 
to the study of groups, but not as good for the 
study of rarer, more massive clusters. 

We identify galaxies using SKID^. Our 60- 
particle galaxy mass resolution limit (5 x IO^^Mq) 
corresponds to « i*/4, based on the equivalent 
number density of observed galaxies in the Sloan 
Digital Sky Survey (Blanton et al. 2001; for 
further discussion of the resolution limit and com- 
pleteness threshold of the simulated galaxy sample 
see Murali et al. 2002). We identify bound sys- 
tems using a spherical-ovcrdcnsity (SO) criterion 
on friends-of- friends halos, as described by Gard- 
ner (2001). In this method, the most bound 
particle within a friends-of-fricnds halo is found, 
and the system extends radially outwards until 
the enclosed average density reaches the Virial 

^ Spline Kernel Interpolative DENM AX, publicly available at 
http : //www-hpcc . astro . Washington . edu/tools/skid . html. 



density threshold (282 times the mean density for 
our cosmology) . SO systems that contain three or 
more galaxies above our resolution limit are iden- 
tified as "groups." We find 128 groups at z = 0, 
spanning a mass range of 10^^-^ — lO^'^'^M©, with 
up to 42 member galaxies. Note that all of these 
groups are bound; we identify groups in three di- 
mensions, so our catalog does not contain chance 
line-of-sight projections. 

To study the effects of numerical resolution, we 
also analyze a smaller volume, higher resolution 
simulation having the same cosmological model; 
we will refer to this as our "high-resolution" sim- 
ulation. This simulation has 2 x 128'^ particles in 
a 22.222/i~^Mpc volume, with 3.5/i~^kpc (equiva- 
lent Plummer, comoving) softening, yielding par- 
ticle masses of mbary = 1-05 x IO^Mq and rridark = 
8.8 X lO^Afg. We find 52 groups in this simula- 
tion at z = 0. By chance, this volume contains 
one anomalously large (10^'''^Mq) group contain- 
ing 217 galaxies above our 60-particle galaxy mass 
limit of 6.3x IO^Mq, but the other 51 groups probe 
the expected mass range down to lO^^'^Af0. 

The velocity dispersion of a group can be cal- 
culated from its dark matter (ctdm) or its member 
galaxies (ugai). We quote 1-D velocity dispersions, 
obtained by dividing the 3-D velocity dispersion 
by \/3. We find that ctdm correlates tightly with 
group mass while dgai shows systematic deviations 
from aDM for smaller groups, as discussed below 
in §5.3. 

To calculate X-ray emission, we use the latest 
plasma code of Raymond & Smith (1977). The in- 
terface to simulations via the TIPSY^ package was 
kindly provided by G. Lewis and C. Murali. We 
calculate particle luminosities in the 0.5 — 2 keV 
bandpass that is most commonly employed to 
study group X-ray properties using ROSAT. Un- 
less otherwise noted, the X-ray temperature is cal- 
culated from the average luminosity- weighted tem- 
perature of group particles; we will discuss some 
issues relating to this in §4.3. 

Applying the usual SPH formalism to obtain 
gas densities for calculating luminosities results 
in a gross overestimation of the X-ray emission. 
The usual SPH algorithm does a poor job of re- 
solving the two-phase interface between the intra- 
group medium and gas in galaxies. As a conse- 

^Publicly available at http://ww-hpcc.astro.washington.edu/tools/tipsy/tipsj 
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quence, the density of a hot intragroup gas par- 
ticle is greatly overestimated when it happens to 
lie very close to a cold, dense clump, producing 
large and unphysical X-ray cmissivitics for some of 
the particles. To correct for this effect, we follow 
Pearce et al. (2000) and explicitly decouple the 
hot and cold phases by recalculating gas densities 
using only hot particles with T > lO^K. Unlike 
Pearce et al., however, we only implement decou- 
pling at the post-processing stage, not during the 
simulation itself. Our results are insensitive to the 
choice of temperature threshold, so long as it is 
below T < IG^K, since virtually all gas particles 
bound in groups that have not cooled into galax- 
ies have temperatures higher than this threshold. 
Croft et al. (2001) performed extensive tests of 
this approach, and found that it is indeed insensi- 
tive to the exact temperature threshold, behaves 
as expected with resolution, and reproduces the 
correct X-ray luminosity in analytic cases. The 
number of particles affected is small, but the to- 
tal correction factors arc large (typically a few to 
ten; see Croft et al. 2001), so our luminosity pre- 
dictions must be read with this caution in mind. 
Additionally, since we do not explicitly decouple 
the hot and cold phases while evolving the sim- 
ulation, we probably overpredict the amount of 
cooled baryons, by an uncertain amount. (The 
cold phase is ultimately "decoupled" from the hot 
phase by conversion into collisionless stars, but 
this conversion is not instantaneous.) Springel & 
Hernquist (2001) suggest that an alternative for- 
mulation of the SPH equations can resolve this 
"decoupling problem" in a self-consistent manner, 
and we will investigate this possibility in future 
simulations. 

3. Testing Self-Similarity 

3.1. Simulated Group Scaling Relations 

We begin by examining whether our simu- 
lated groups follow the self-similar scaling rela- 
tions expected for systems with negligible non- 
gravitational heating. To do this, we compute 
X-ray luminosities using an intragroup metallic- 
ity of zero, we sum the luminosity of all gas in 
the group out to its Virial radius, and we take 
the group's velocity dispersion to be that of the 
dark matter. If group formation were well de- 
scribed by non-radiative physics (a.k.a. "adia- 



batic" physics, though gravitational shock heat- 
ing is also included), then this approach should 
closely reproduce the self-similar scaling relations 
(see, e.g., Owen et al. 1998; Muanwong et al. 
2001). 

Figure 1, upper left, shows the Lx — ctdm re- 
lation for our groups. The plot symbols indicate 
the number of galaxies in each group. A break is 
clearly evident, occurring at ctdm ~ 180 km/s; the 
fit above the break is extended to lower ctdm as 
the dashed line, and it is clearly discrepant with 
these smaller systems. The dotted line shows the 
self-similar scaling relation, Lx cx tTpj^, normal- 
ized to our largest groups. The break value is 
determined by requiring continuity between the 
fits above and below the break. Best fit power 
law relations above and below this break value are 
listed in the lower right. Even above the break, the 
Lx ~ f DM relation is steeper than predicted by the 
self-similar model, with Lx oc (J^l^^'^- Below the 
break the relation is very steep, Lx oc cr^l^^'^- 
The fits are unweighted, to follow the procedure 
most commonly used by observers. We list our 
scaling relation fits in Table 1, including the scal- 
ing relations observed by MZ98 and HPOO. We 
defer a comparison with observations until §5.4, 
where we include various observational effects in 
our analysis; here we simply note that the observa- 
tions show a similar qualitative trend in which lu- 
minosities drop faster with decreasing group mass 
than the self-similar model predicts. 

Figure 1, upper right, shows the Lx — Tx re- 
lation. Here again a break is evident, though less 
obviously so, at T « 0.7 keV. Above the break, 
the relation is Lx cc tJ'^^"'^, clearly discrepant 
with self-similarity (Lx oc T^, dotted line), while 
below the break it steepens to Lx oc T^'^^°'^. A 
single power-law fit yields a slope of 3.4±0.1. Here 
the break value is somewhat more ambiguous, but 
even with the break placed at T = 2 keV (in which 
case the fits are clearly discontinuous), the slope 
above the break is still 2.6, and inconsistent with 
self-similarity. Though a double power law is an 
adequate fit, there is an overall trend for a contin- 
uous steepening of the relation to lower tempera- 
tures. 

Our luminosity scaling relations demonstrate 
two points. The first is that a break occurs at 
all, when the simulation incorporates no physics 
"tuned" to pick out this mass or temperature 
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Fig. 1. — Scaling relations of the simulated galaxy 
groups, assuming zero metallicity. The upper left 
panel shows Lx — cdm, upper right Lx — Tx, and 
lower left Tx — ctdm- The symbols indicate the 
number of galaxies in each group. The best-fit re- 
lations are listed in the lower right. The solid line 
shows the best fit, having a break at the value indi- 
cated in the legend (except in the lower left panel, 
where a single power law fit is acceptable). The 
dashed line is a continuation of the slope above 
the break to lower values of ctdm or Tx ■ The dot- 
ted lines show predictions from self-similar scaling, 
normalized to our largest groups. These scaling 
relations may be compared to, e.g., Figures 4-6 of 
Mulchaey (2000); we will present a more detailed 
comparison to observations in §5.4. In no regime 
do our simulated luminosity scaling relations fol- 
low those predicted by self-similar models and by 
adiabatic numerical simulations. 



scale. The second is that even above the break, 
the slope shows a significant departure from self- 
similarity. In no regime do our groups follow 
the expected self-similar relations, suggesting that 
self-similarity is an inappropriate model for relat- 
ing clusters and groups. This strong departure 
from self-similarity is our most important result; 
the rest of the paper will be devoted to under- 
standing its origin and comparing our predictions 
to observed group properties. 

The lower left panel of Figure 1 shows the 
Tx — ctdm relation. In contrast to the luminos- 
ity relations, this relation shows no break, and the 
scaling, T oc al^^'^° °^, is only slightly off from the 
self-similar prediction {Tx oc tr^^, dotted line). 
This agreement suggests that the X-ray tempera- 
ture and dark matter velocity dispersion are both 
determined mainly by the gravitational potential, 
as expected for systems that are Virialized and in 
hydrostatic equilibrium. The slight excess heating 
that makes the relation shallower than Tx oc a-^y^ 
may be coming from shock heating of gas on fil- 
aments before accretion into the group potential, 
which can heat gas to ~ 10^ K (Dave et al. 2001). 
Supernova feedback as implemented in this simu- 
lation is not expected to add any significant heat 
to the intergalactic gas. In any case, the amount 
of heating indicated in this plot is far from suffi- 
cient to explain the departures from self-similarity 
in the Lx — Tx relation. Moreover, since there is 
no evidence of a break in Tx — ctdm, it is clear 
that some other physical process is at work in the 
luminosity scaling relations. 

3.2. Physical Origin of the Simulated Scal- 
ing Relations 

Why does our simulation yield lower-than- 

expected luminosities as one goes to smaller sys- 
tems? The origin of this effect does not appear to 
be some heat source that adds pressure support 
and lowers the central density at a given Virial 
mass, nor does it appear that the low mass sys- 
tems have lower luminosities because of substan- 
tially sub- Virial temperatures. Either of these 
scenarios would require a feature in the T^ — a 
relation close to the break at <t « 180 km/s in the 
Lx — cr relation, and none is seen. Alternative ex- 
planations are that the X-ray emitting gas density 
is lowered by removal of gas from this phase, or 
that the density structure of the gas is changing 
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Table 1 

Simulated Group Scaling Relations. 





SB'^ 


cr 


r d log Lx 1 


rdlogix 1 


rdlogix 1 


rdlogix 1 


rdiogTxi 


L dlogo- J >180 


L dlogo- J<180 


LdlogTx J >0.7 


LdlogTx J <0.7 


L d log a J 





none 




4.7 ±0.3 


9.7 ±0.6 


2.9±0.1 


4.1 ±0.5 


1.86 ±0.05 


eq. 6 


none 




4.3 ±0.2 


9.5 ±0.6 


2.5 ± 0.2 


4.6 ±0.3 


1.69 ±0.05 




none 


cdm 


3.8 ±0.2 


7.7 ±0.5 


2.2 ±0.2 


4.0 ±0.3 


1.77 ±0.05 





eq. 1 




5.0 ±0.3 


11.9±0.9 


3.4±0.2 


4.8 ±0.4 


1.88 ±0.06 





none 


0-gal 


3.9 ±0.4 


3.9 ±0.9 


2.9 ± 0.1 


4.1 ±0.5 


1.14 ±0.07 


eq. 6 


eq. 1 




4.4 ±0.2 


11.7±0.9 


2.6 ±0.1 


4.8 ±0.4 


1.75 ±0.05 


MZ98 






4.3 ±0.4 




2.8±0.1 




2.2 ±0.9 


HPOO 






4.5 ± 1.1 




4.9 ±0.8 




0.9 ±0.2 



'^Assumed metallicity. 
^Surface brightness correction. 



systematically with group mass. Both of these ef- 
fects are, in fact, occurring in our simulation, and 
together they account for the departure from the 
self-similar scaling relations. 

The top panel of Figure 2 shows the fraction 
/hot = Miot/Mbaryon of cach group's total baryons 
that are in the hot, X-ray emitting phase. There 
is a clear trend that smaher systems have less 
X-ray emitting gas. Our largest systems have 
slightly more than half of their baryons in a hot 
phase, while our smallest systems have less than 
20%. The trend steepens slightly to smaller sys- 
tems. The steepening is visually exacerbated when 
one plots /hot^ on a logarithmic plot, hence the 
appearance of a "break" in the Lx — Tx and 
Lx — cdm relations. 

The data points in Figure 2 are binned from 
a compilation of observations of 30 X-ray groups 
and clusters with temperatures ranging from 0.8- 
4.5 keV. The lowest 17 are from Mulchaey et al. 
(1996), the next five are from Hwang ct al. (1999), 
and the highest eight are T < 4.5 keV Abell clus- 
ters from Cirimele, Nesci, & Trevese (1997). The 
error bars shown indicate the error in the mean 
computed from the scatter among the data points 
in each bin. In each case, we take the ratio of the 
gaseous to gaseous±stellar mass estimates. The 
observations do show a trend of decreasing hot 
gas fraction with decreasing group mass. Unfor- 
tunately, current data do not extend to the mass 



range of the lowest mass groups in our simulation, 
where the reduction in hot gas is largest. 

In the top panel of Figure 2, there appears to 
be remarkably good agreement between the sim- 
ulated and observed hot gas fractions. Unfortu- 
nately, this agreement is somewhat misleading. As 
Balogh et al. (2001) have argued. X-ray observa- 
tions of groups do not extend nearly as far out, 
relative to the virial radius, as they do in clusters, 
and since the hot gas is distributed more diffusely 
than the cold gas, the observed hot gas fraction in 
smaller groups is biased low. While it is difficult 
to correct the observations, it is possible to model 
this bias in our simulated groups to enable a fairer 
comparison. To do so, wc recalculate the baryon 
fractions only within a radius given by 

where Rx is the extent of the X-ray emission. 
This relation is obtained by fitting the data in 
Figure 3 of Mulchaey (2000) [for which we ob- 
tain Rx/Rvir = (logr + 0.7)/1.7] and using the 
Tx — relation from Figure 1 . Wc take Ryi^ to 
be the maximum extent of the group. 

The result of applying equation (1) is shown in 
the second panel of Figure 2. The predicted hot 
fraction is indeed smaller than in the top panel, in- 
creasingly so for smaller systems. Thus, in a fairer 
comparison with data, our simulation appears to 
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Fig. 2.— Top panel: Ratio of hot (T > lO^K) 
to total baryons as a function of group veloc- 
ity dispersion. The top axis legend indicates the 
approximate group temperature, obtained by us- 
ing the Tx — tJDM scaling relation found in Fig- 
ure 1. The data points show binned observa- 
tional values as described in the text. Second 
panel: Same as top panel, but results have been 
computed only out to an "observable" fraction 
of the Virial radius, as given by equation (1). 
Third panel: Same as second panel, but based 
on the high-resolution, smaller volume simulation. 
Fourth panel: "Clumping factor" (eq. 2) of hot gas 
in each group, from the large- volume simulation. 



"ovcrcool" baryons in groups. As wc have dis- 
cussed elsewhere (Dave et al. 2001; Murali et al. 
2002), the simulation also appears to "overcool" 
baryons globally the fraction of baryons con- 
verted to stars is 24%, which is high compared to 
estimates from the observed luminosity function 
(Blanton et al. 2001; Cole ct al. 2001), though it 
is consistent with estimates from the extragalactic 
background light (Bernstein, Freedman, & Madore 
2001). 

Pearcc ct al. (2000) and Balogh et al. (2001) 
have argued that a higher resolution simulation 
should produce substantially more cooled gas. We 
can test this expectation directly using our high- 
resolution simulation. The third panel of Fig- 
ure 2 shows the hot fraction for the 52 groups 
identified in this simulation, including the surface 
brightness cut described above. Interestingly, the 
trend is opposite to that argued previously: there 
is actually somewhat more hot gas in the higher 
resolution simulation. We believe that this in- 
crease in hot gas fraction with increasing numeri- 
cal resolution, i.e. a decrease in galaxy masses, is 
caused by the difficulty in resolving the interface 
of two such distinct phases. This clauses cooling 
rates in the immediate vicinity of galaxies to be 
overestimated (sec §2), and affects a smaller frac- 
tion of particles in the higher resolution simula- 
tion. The hot gas fractions in our high-resolution 
simulation, including surface brightness bias, are 
in good agreement with the observations, though 
the range of overlap is still limited, and we do 
not have a still higher resolution simulation with 
which to demonstrate convergence of the numeri- 
cal prediction. The global cold fraction in the high- 
resolution simulation is slightly higher than that 
of the large-volume simulation, with 26% of the 
baryons converted to stars, since the simulation 
resolves cooling and star formation further down 
the mass function. An estimate of the total stel- 
lar mass fraction, obtained by combining together 
several simulations of varying volume and resolu- 
tion is presented in Fardal et al. (2002 in prepa- 
ration). 

The simple uniform removal of hot gas is not 
sufficient to fully explain the reduced luminosi- 
ties of lower mass groups. Consider groups with 
Tx « 3 keV (logT « 0.5), having hot gas frac- 
tions ~ 0.5, and those with Tx ~ 0.3 kcV, having 
hot gas fractions ~ 0.2. Relative to a model with 
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constant gas fraction, the low mass groups in our 
simulation have hot gas densities lower by a factor 
of 2.5. Since Lx oc n^ot: this reduction should de- 
press the luminosities of the low mass groups by 
a factor of 6.25. However, the gap between our 
best-fit relation and the self-similar scaling (solid 
and dotted linos of Figure 1, at logTx = —0.5) 
is approximately a factor of 25. This shows that 
the uniform removal of gas accounts for somewhat 
more than half (logarithmically) of the reduction 
in luminosity. 

The remaining factor arises because hot gas is 
less clumped in the smaller systems. The bottom 
panel of Figure 2 shows the clumping factor 

(Phot) A'hot/Ohot 

for each of our systems. The (..) in the first 
equality represent volume averages, and the sec- 
ond equality expresses this ratio as a sum over 
(equal mass) SPH particles, where pi is the SPH- 
estimated density of particle i, TVhot is the number 
of hot gas particles in the group, and phot is the 
mean density of hot gas, namely the total hot gas 
mass of the group divided by the volume out to the 
Virial radius. If the hot gas in these groups were 
isothermal, then this clumping factor would be di- 
rectly proportional to the luminosity (cf. Evrard et 
al. 1996), since the volume cmissivity is propor- 
tional to p^Qj throughout the group. If Chot were 
the same in high and low mass groups, then their 
relative luminosities would scale in proportion to 
/hot (plus temperature factors, but the group tem- 
peratures do follow the self-similar scaling to a 
good approximation). The term "clumping fac- 
tor" should be applied with some caution, since 
the globally averaged quantity defined in equa- 
tion (2) could be different in two groups that have 
smooth gas distributions but different radial den- 
sity profiles Phot('')- 

Figure 2 shows that Chot drops by a factor of 
~ 3 from 3 keV systems to 0.3 keV systems (ac- 
cording to the best-fit relation shown as the solid 
line), though the scatter is large. This indicates 
that clumping is responsible for most of the re- 
maining drop in luminosity, after accounting for 
changes in /hot- The physical reason for the drop 
in Chot is not obvious. One possibility is that 
larger systems are younger and hence less dynami- 
cally relaxed, and therefore contain more substruc- 



ture. Another possibility is suggested by the an- 
alytic model of Bryan (2000), who also predicts 
a luminosity drop of a similar factor due solely to 
reduced "clumping" in smaller systems (compare 
his dot-dashed and solid lines in his Figure 3). In 
his case, this difference arises because cooling more 
efficiently removes hot gas in the central regions of 
smaller systems. As we will show in §4, we do not 
see a drop in the mean density in the inner regions 
of smaller systems, however the differences in Chot 
do arise mainly in the inner regions. Regardless of 
cause, it is the combination of lowered subclump- 
ing and an overall drop in the hot fraction that is 
responsible for the lowered luminosities in smaller 
groups. 

Figure 3 shows the relative impact of these two 
effects on the Lx — Tx relation for groups in our 
simulation. We focus on Lx — Tx here because 
its deviation from self-similarity is more evident 
than that of Lx — cdm, at least above the break 
(cf. Figure 1). The solid points (at the bottom 
of the vertical lines) represent the groups repro- 
duced from Figure 1. The open circles represent 
the group luminosities "corrected" to a common 
gas fraction and clumping factor, 

-^X.corrected = -^x(0.5//hot)^(C3 kev/Chot), (3) 

where /hot is the fraction of hot baryons in each 
group, C3 keV is the typical clumping factor for 
3 keV groups, and Chot is the clumping factor for 
the specific group (cf. Figure 2). If all groups 
had a hot fraction of 50% and a clumping factor 
of « 20, this correction would have no effect. In- 
stead, there is an increasingly larger correction for 
smaller systems. The solid portion of each verti- 
cal line represents the correction factor due to the 
mean hot fraction, and the dotted portion repre- 
sents the remaining factor due to clumping. Fig- 
ure 3 shows that the two effects are of compa- 
rable magnitude, and that the combination of the 
two brings the group luminosities into much better 
agreement with the self-similar relation (diagonal 
dotted line), though there is still some tendency 
for the lower mass groups to be under luminous. 
The best-fit relation for the "corrected" luminosi- 
ties (which are not, of course, the physically pre- 
dicted luminosities) has a slope of 2.4±0.2. The re- 
maining discrepancy with the self-similar Lx — Tx 
relation, and some of the remaining scatter, proba- 
bly arise from the temperatures themselves, since 



9 




log (keV) 

Fig. 3. — Lx — Tx relation showing the relative 
importance of variations in mean hot gas density 
and clumping factor among our groups. The solid 
points are reproduced from Figure 1 (but not la- 
beled by the number of member galaxies), while 
the open circles include correction factors to force 
all groups to have the same mean hot gas fraction 
and clumping factor (see eq. 3). Self-similar scal- 
ing is shown as the dotted line, and the best fit to 
the open circles is shown as the solid line, having 
a slope of 2.4 ± 0.2. The two portions of the ver- 
tical line segments represent the separate correc- 
tions due to hot gas fraction (solid) and clumping 
(dotted). 



smaller groups tend to be hotter than the self- 
similar Tx — cr relation predicts (see Figure 1). 

A natural concern is that the clumping factor of 
our smaller groups is artificially suppressed by res- 
olution effects, since our smallest systems contain 
as few as ~ 100 hot particles. Figure 4 shows a 
comparison of Lx — ctdm from our high-resolution 
simulation (open circles) to the data shown in 
Figure 1 from our large-volume simulation (solid 
points). We focus on Lx — (Jdm because ctdm is 
a more direct indicator of group mass than Tx- 
While the mass range probed by the two simula- 
tions is necessarily different, in the region of over- 
lap they agree remarkably well. The best fit rela- 
tion to the combined data is shown in the lower 
right, and is in good agreement with the scalings 
from Figure 1. If resolution were a major fac- 
tor in the predicted departures from self-similar 
scaling, we would expect low mass groups to be 
systematically more luminous in the higher reso- 
lution simulation, and we would expect the break 
at (TDM ~180 km/s to shift down by a factor of 
two to Ctdm ~90 km/s, corresponding to the factor 
of eight increase in mass resolution. While there 
is only one point well above the break from the 
high resolution simulation, it sits close to the mean 
Lx — cdm relation defined by the lower resolution, 
large-volume simulation. The good agreement in 
Figure 4 demonstrates that the reduced luminos- 
ity of low mass groups in our simulations is not an 
artifact of finite resolution but arises instead from 
the physical processes that the simulations incor- 
porate. This test also indicates that our two-phase 
decoupling method for calculating luminosities is 
robust against modest changes in resolution. 

In summary, the breaks in the Lx — ctdm and 
Lx — Tx relations are driven by the increased ef- 
ficiency of radiative cooling in lower mass groups, 
which affects both the fraction and the density 
structure of the hot gas. The fraction of T > 10'' K 
gas for our simulated groups drops from 50% at 
a w 500 km/s to 20% at ct w 100 km/s, which 
explains the majority of the decline in luminos- 
ity relative to the self-similar prediction, while the 
remaining reduction is due to smaller systems hav- 
ing lower clumping factors Chot ■ This trend in hot 
fraction is consistent with observations, though 
the hot gas fractions of our large-volume simula- 
tion appear somewhat too low once surface bright- 
ness biases are taken into account. The hot gas 
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Fig. 4. — Resolution test of the Lx — ctdm rela- 
tion, showing the groups from the (50/i~^Mpc)^ 
simulation (solid points, as in Figure 1 but not 
labeled by the number of member galaxies), and 
the higher resolution, (22.22/i~^Mpc)'^ simulation 
(open circles). The agreement in the region of 
overlap is very good, indicating that our results 
for the Lx — ctdm relation are not significantly 
affected by numerical resolution. The scaling re- 
lations of the combined sample are shown as the 
solid line, and listed in the lower right. The dot- 
ted line shows self-similar scaling normalized to 
the largest systems. 



fractions in our high-resolution simulation agree 
very well with observationally inferred values, and 
the Lx — CTDM relations of the two simulations 
agree very well in their range of overlap. 

The trend of a declining hot fraction in the 
group mass regime seems to be a generic feature 
of CDM-based galaxy formation, arising in hydro- 
dynamic simulations of galaxy formation (Kay et 
al. 2000; Blanton et al. 2000; Pearce et al. 
2001; this paper) and semi-analytic models (Bower 
et al. 2001). At the present epoch, the cool- 
ing time in the intragroup medium is long, but 
during the early stages of their assembly, a larger 
fraction of baryons were able to cool in smaller 
systems. This rather straightforward physical pro- 
cess just happens to become important in the mass 
regime of groups, and it results in a breaking of 
self-similarity in the observed scaling relations. 

4. Profiles 

Radial profiles of physical quantities such as 
temperature, surface brightness, and electron den- 
sity can provide insights into galaxy formation 
processes within groups. In this section we present 
our simulation results for these profiles, fit beta 
models to the surface brightness profiles, and com- 
pare with observations, including the entropy- 
temperature relation derived in PCN99. 

4.1. Average Profiles 

Figure 5 shows the average surface bright- 
ness, temperature, electron density, entropy, and 
clumping factor profiles for our simulated groups, 
divided into four mass bins with 32 groups per 
bin. Profiles are scaled to the group Virial radius, 
then averaged. The surface brightness profiles 
(top panel) are shown scaled to the group central 
surface brightness. They are generally similar in 
form down to our smallest groups. This result 
contrasts with, e.g.. Figure 1 of PCN99, where 
the surface brightness profile is significantly shal- 
lower for smaller systems. However, MZ98 find 
the shapes of surface brightness profiles for poor 
groups to be quite similar to those of clusters, 
consistent with our simulated profiles. We will ex- 
amine these issues in more detail when we discuss 
beta models in the next section. 

Our temperature profiles (second panel) show 
a slight drop from the inner to outer parts. We 
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Fig. 5. — Radial profiles of surface brightness 
(top panel), temperature (second), electron den- 
sity (third), entropy (fourth), and clumping fac- 
tor (bottom) for the simulated groups, divided by 
mass into four equal number subsamples (solid, 
dotted, short-dashed, and long-dashed lines, from 
highest to lowest mass). The surface brightness 
and temperature profiles are projected, averaged 
in cylindrical shells, while the other profiles are 
3-d, averaged in radial shells. 



do not produce a cool central region often seen in 
groups (e.g., HPOO) or cooling flow clusters (e.g., 
Allen, Schmidt, fc Fabian 2001); we will return 
to this issue in §4.3. However, beyond the in- 
nermost regions, both group and cluster temper- 
ature profiles are observed to have a slight drop 
with radius (MZ98; HPOO; dc Grandi & Molendi 
2002). In clusters, where one is able to probe al- 
most out to the Virial radius, BeppoSAX obser- 
vations by dc Grandi & Molendi (2002) indicate 
an isothermal core (for non-cooling flow clusters) 
out to 0.2i?vir, and then a 30-40% drop out to 
0.5i?vir, which they claim is a steeper drop than 
predicted by numerical simulations. From 0.2i?vir 
to 0.5i?vir, our Smallest groups show a 20% drop in 
temperature, while our larger groups show a 40% 
drop, consistent with their observations. Out to 
the Virial radius, the drop is around a factor of 
two. Thus, with minor discrepancies, our simu- 
lations generally reproduce observed temperature 
profiles outside of 0.2i?vir, predicting temperature 
profiles that are falling slowly towards large radius. 
We have presented projected temperature profiles 
in Figure 5 to allow more straightforward compar- 
ison to observations, but our conclusions from the 
3-D profiles would not be significantly different. 

The electron density profiles (middle panel) for 
the three subsamples would lie on top of each 
other in a purely self-similar model. The fact that 
smaller groups have lower rig is a signature that 
cooling has operated more efficiently in these sys- 
tems, as we showed in Figure 2. Figure 5 shows 
that cooling preferentially removes hot gas at in- 
termediate radii. It is here that the cooling time 
is comparable to a Hubble time, and hence there 
is the most dramatic difference between large and 
small systems. In the central region, the cooling 
time has been short compared to a Hubble time 
over much of the group's life, while in the outskirts 
it is much longer than a Hubble time; thus in both 
regimes the electron density is similar across all 
groups. The shapes of our simulated density pro- 
files arc markedly different from those predicted 
in a pre- heating scenario, where the electron den- 
sity profile would show an increasingly large core 
in smaller systems. Unfortunately there are no 
direct observations of the electron density profile, 
only estimates from beta-model fitting of surface 
brightness profiles (PCN99). The combination of 
rapidly improving Sunyaev-Zel'dovich and X-ray 



12 



observations of groups may yield electron density 
profiles in the future. 

The entropy profile (second from bottom) is 
computed by taking the mass-averaged value of 
T/n^J^ of all particles within each bin. The en- 
tropy profile declines smoothly from the outer edge 
into the center, in agreement with simulations 
by Thomas et al. (2002) and analytic expecta- 
tions based on convective equilibrium (Voit et al. 
2002), though they differ from those in the ana- 
lytic model of Babul et al. (2002). We will discuss 
entropy further in §4.4. 

The clumping factor profile (bottom panel) 
shows Pi/NhotPhot (cf. eq. 2, but now the sum 
is over particles in each radial bin, and A^hot and 
Phot are the particle number and mean density of 
hot gas in the bin). The largest systems show 
a noisy profile in the outer parts, presumably 
because these objects are dynamically younger 
and have more subclumps due to recent merg- 
ing. However, the systematic mass dependence of 
these profiles arises mainly at small radii, where 
larger systems show significantly more clumping. 
It is this inner region that drives the trend seen 
in the bottom panel of Figure 2, and consequently 
(a portion of) the luminosity drop in Figure 3. 
In part this trend arises because the largest sys- 
tems often contain a single dominant galaxy in the 
center, whereas in smaller systems there can be 
significant offsets between the group center and 
the largest galaxy's position. The lower Chot val- 
ues at small radii in smaller groups do not appear 
to be resolution artifacts, since the luminosities 
predicted by the high-resolution and large- volume 
simulations agree (Figure 4) and the electron den- 
sity profiles do not fiatten in the centers of small 
groups, as would be expected if resolution effects 
were important. 

In summary, profiles may have greater power 
to discriminate between models of group and clus- 
ter formation than scaling relations alone, since 
the latter can be adequately reproduced in both 
cooling and pre- heating models. In particular, the 
surface brightness profile appears to be a key di- 
agnostic, with initial indications being that our 
simulations may be in disagreement with observa- 
tions. We now examine this issue in more detail. 



4.2. Beta Models 

The canonical form for representing an X-ray 
surface brightness profile is a "beta model". In 
this model, the parameter (3 represents the ratio 
of the specific cniergy in dark matter to that of the 
hot gas, namely 



where /x is the molecular weight and rup is the pro- 
ton mass. For an isothermal sphere in hydrostatic 
equilibrium with a King model density profile, the 
resulting surface brightness profile is 

5(r) = 5(0)(l+r/reore)"'''+°-', (5) 

where S{r) is the azimuthally averaged surface 
brightness at projected radius r, and rcore is the 
core radius. The electron density distribution also 
follows a similar form, with an exponent of —3/3/2. 

Observationally, /3 may be estimated either via 
equation (4), inferring the gas temperature from 
the X-ray spectrum and obtaining the velocity dis- 
persion or mass independently, or via equation (5), 
by fitting surface brightness profiles. These spec- 
troscopic and imaging estimates of /3 arc usually 
referred to as /Jgpec and (3^, respectively. Unfortu- 
nately, the two estimates have historically yielded 
discrepant answers even in well-studied clusters, 
with /3fit « 2/3 (Mohr, Mathiesen & Evrard 1999) 
and /?spec ~ 1 (see discussion in Mulchaey 2000). 
Using hydrodynamic simulations, Navarro, Frcnk, 
& White (1995) find that fitting surface bright- 
ness profiles yields pat ~ 0.8 even when /?spec ~ 1 
because the dark matter mass profiles deviate sys- 
tematically from a King model. Furthermore, 
they demonstrate that fitting to a smaller frac- 
tion of the Virial radius (as often done in obser- 
vational analyses, especially of groups) biases Pat 
even lower. These biases were invoked to rec- 
oncile observations of /3spcc and Pat- However, 
recent Chandra observations of six large clusters 
with lensing masses suggest Pspec ~ 0.7 (Allen, 
Schmidt, & Fabian 2001), so there may be no 
discrepancy to be explained after all. 

Determinations of Pat for groups are even more 
unsettled, primarily because only a small por- 
tion of the surface brightness profile is visible in 
these faint systems, and also because the contri- 
bution from a central galaxy, which is excluded 
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to obtain the intragroup gas properties, is typ- 
ically more important than in clusters. HPOO 
find a median /?fit = 0.46 from 24 groups with 
Tx ~ 0.4 — 1.6 kcV. This implies considerably 
flatter surface brightness profiles than in clusters, 
a fact that is the basis for the "entropy floor" of 
PCN99. However, MZ98 find a median /3fit » 0.8 
for 9 groups, which is higher than typical observed 
values for clusters. Both MZ98 and HPOO use two- 
component beta models that are statistically war- 
ranted by the data, and in fact they consider many 
of the same groups. The fits of HPOO for groups in 
common with MZ98 yield much smaller core radii 
and smaller /3fit, suggesting a strong degeneracy 
between rcore and /?fit (Mulchaey 2000). HPOO 
claim to sec a trend for lower /3fit in smaller sys- 
tems, but this could be an artifact of the bias men- 
tioned above, where fits to a smaller fraction of 
the Virial radius yield lower /9fit values, since ob- 
servations of their smallest groups extend to only 
~ 10 — 20% of the Virial radius. Thus, in our view, 
the case for flatter surface brightness profiles in 
smaller groups has not been clearly proven. 

Figure 6 shows a histogram of /3fit values from 
two-component beta model fits to our simulated 
group surface brightness profiles. The four line 
types correspond to our equal number subsam- 
ples. All groups with /3fit > 1.6 are stacked in the 
rightmost bin. We allow the two components to 
simultaneously vary freely in our Marquardt mini- 
mization routine, then take the (3 value associated 
with the component having the larger Tcoro, since 
it presumably represents the intragroup compo- 
nent^. The median values of /Jgt and rcore are 
0.66 and 0.2i?vir, respectively. The second com- 
ponent (not shown) has median (3f^t « 0.5 and 

^core ^ 0.03-??vir. 

Our median = 0.66 is intermediate between 
that of HPOO and MZ98. We find only a slight 
trend with size: our subset of largest groups has 
a median /Sgt = 0.68, while our subset of small- 
est groups has median /3fit = 0.62. We find that 
for many groups, the two-component beta model 
is not a good fit even though sensible values of 
/3fit are returned. Since groups are not strictly 



^In the case of groups, where the central galaxy may make 
a significant contribution to the gravitational potential, we 
are not convinced that it is sensible to fit the "galaxy" as 
a physically distinct component, but here we follow the 
standard practice of observational analyses. 
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Fig. 6. — Top panel: Histogram of (3^t for the 
intragroup surface brightness component of two- 
component beta models fit to our groups. Median 
/Jfit value is indicated in upper right. Middle panel: 
Histogram of /3spec for luminosity-weighted group 
temperatures. Bottom panel: Histogram of /3spec 
for mass-weighted average group temperature. In 
each panel, the four histograms correspond to 
four mass bins containing nearly equal numbers 
of groups, with solid, dotted, short-dashed, and 
long- dashed lines in decreasing order of mass. 
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isothermal, do not follow King model profiles, and 
are typically not spherically symmetric (especially 
smaller systems), perhaps this result is not surpris- 
ing. It may, however, explain why the observed 
fits are so sensitive to the details of the fitting 
procedures. Further evidence that beta models 
arc inappropriate comes from a comparison of the 
surface brightness and electron density profiles in 
Figure 5; though the former have a similar shape 
across all systems, the latter do not. 

In principle, surface brightness profiles offer an 
excellent diagnostic of intragroup gas physics. Un- 
fortunately, current observational determinations 
of Pat appear too uncertain to validate or in- 
validate any models. This situation is likely to 
change with XMM and Chandra, given their in- 
crease in sensitivity and spatial resolution. If ob- 
servations of flatter surface brightness profiles in 
smaller groups persist (i.e., if Pat correlates with 
temperature), they would pose a serious challenge 
for our simulation predictions. Similar values of 
/3fit in groups and clusters, on the other hand, 
would strongly favor the scenario presented here, 
in which cooling rather than pre- heating is the pri- 
mary cause of deviations from self-similar scaling. 
Directly measured profiles will be more informa- 
tive than model fits, since beta models implicitly 
incorporate a number of assumptions that may not 
be true in detail. 

4.3. Group Temperatures cind Pspec 

The middle and bottom panels of Figure 6 show 
histograms of Pspec from our simulated groups. 
The middle panel shows Pspec when the usual 
luminosity-weighted temperature is used, and the 
bottom panel shows the case when the mass- 
weighted temperature is used; in both cases we 
have used ctdm for the velocity dispersion. In- 
terestingly, the luminosity-weighted case shows 
significant deviations from the expected /? = 1, 
with a median value of /Jgpcc = 0.47. Thus, 
the luminosity-weighted temperatures are roughly 
twice the groups' Virial temperatures. Conversely, 
the mass-weighted temperatures underestimate 
the Virial temperature by ~ 20%. There is no 
systematic trend of Pspec with group mass in ei- 
ther case. 

Adiabatic simulations yield /3spec ~ 1 (e.g. 
Navarro, Frenk, & White 1995). As emphasized 
by Allen, Schmidt, & Fabian (2001), this value 



appears high compared to the latest Chandra ob- 
servations of clusters, where /3spec ~ 0.7. Hints 
of this discrepancy were already present in ear- 
lier ROSAT data (see Thomas et al. 2002, for 
a summary). In the simulations of Thomas et 
al. (2002), either cooling or pre- heating pro- 
duces increased temperatures compared to the adi- 
abatic case, yielding better agreement with the 
data. With pre-heating, the reason for the in- 
creased temperature is obvious, but in the case 
of cooling, this nonintuitive result arises because 
gas is cooled out of the center of the cluster, and 
the lowered pressure support draws in hotter gas 
from the outskirts (Pearce et al. 2000). The clus- 
ter simulation with cooling by Lewis et al. (2000) 
produced a similar result. In both cases, however, 
the derived value of /Jgpec was consistent with the 
latest observations, while our groups still appear 
too hot. 

We conclude that our predicted luminosity- 
weighted temperatures are somewhat at odds with 
observations. We investigated whether the proce- 
dure of fitting Raymond-Smith models to coarsely 
binned X-ray spectra yielded some systematic 
temperature bias, and found that when groups 
have a temperature gradient towards the center 
(as in our larger systems), the fitted temperature 
underestimates the true luminosity-weighted tem- 
perature. However, the magnitude of this effect is 
quite small (~ 10%), and in any case it cannot ex- 
plain the discrepancy in our smaller groups, which 
are nearly isothermal. 

Instead, we believe the source of the discrep- 
ancy lies in the temperature profiles of our sim- 
ulated groups (cf. Figure 5). In the sample of 
HPOO, about half the groups show a clear drop in 
temperature, by as much as a factor of two, within 
^ 10—20% of the Virial radius, and the majority of 
remaining groups are at least consistent with hav- 
ing a cooler central region. Similar results are seen 
in large clusters, though the drop in temperature 
is seen out to a smaller fraction of the Virial radius 
(Allen, Schmidt, & Fabian 2001). Our simulated 
groups, however, show no such drop, and since 
the emission is dominated by the central region, 
this leads to a much higher luminosity-weighted 
temperature. If we alter the temperatures of sim- 
ulated group particles "by hand" so that they are 
lower in the center by an amount consistent with 
observations, then our median Pspec increases to 
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roughly 0.8, in agreement with observed values. 
(Specifically, we apply a linear temperature gradi- 
ent such that the central temperature is half that 
at 0.2i?vir.) 

It is not clear why our simulated groups do not 
have a cool central region, while observed groups 
do. If the physical process(es) that produce such 
a region are analogous to those in clusters, then 
the insights gained from recent cluster studies may 
provide clues. For instance, recent observations of 
cooling flow clusters suggest local isothermality in 
the central region (Bochringer et al. 2001a), with 
little emission from the cold component (Peter- 
son et al. 2001), in conflict with the standard 
cooling flow model (Fabian 1994). These findings 
suggest that some unidentified heat source must 
counteract the cooling process, and serve to main- 
tain the central gas in a quasi-static, slightly cooler 
state. This heat soiirce may be AGN in the central 
galaxy (Boehringer et al. 2001b), heat conduc- 
tion from the outer region (Narayan & Medvedev 
2001), outflow induced shocks that can be seen in 
recent Chandra images, or some as yet undiscov- 
ered process. For instance, AGN are observed to 
produce "bubbles" in the intracluster medium, the 
walls of which contain denser, cooler gas (Nulsen 
et al. 2002); any emission weighted measure of 
temperature will be biased low by this gas since 
the emission depends on the density squared. If 
this or similar processes are operating in the cen- 
ters of groups (at a proportionally reduced am- 
plitude that would make them difficult to detect 
directly), then our simulations would be missing 
the input physics required to properly model the 
central regions in these systems. Alternatively, we 
note that the outer radius of the cool region is 
roughly the radius at which the metal-linc; cooling 
rate becomes smaller than a Hubble time (while 
with only the primordial cooling assumed in our 
simulation, the cooling time is still much longer). 
Hence it is possible that including metal-line cool- 
ing in our simulation would result in a cooler cen- 
ter, particularly if there is a metallicity gradient 
towards the center. However, as noted earlier, 
increased cooling does not necessarily reduce gas 
temperatures, since reduced central pressures al- 
low hotter gas to flow in from the surrounding re- 
gions and experience adiabatic heating along the 
way. 



4.4. Entropy 

Entropy is a powerful diagnostic of intragroup 
gas, with a closer connection than temperature or 
density alone to the processes that determine its 
physical state (Bower 1997; Voit et al. 2002). 
The fact that groups appear to show excess en- 
tropy relative to that predicted purely from grav- 
itational shock heating (PCN99) has been used to 
argue for the presence of non-gravitational heating 
processes, since radiative cooling reduces rather 
than increases gas entropy. However, a more com- 
plete analysis suggests that cooling can effectively 
raise the entropy inferred from X-ray data because 
the gas that cools onto galaxies is no longer observ- 
able in the X-rays at all; more efficient cooling in 
smaller systems can therefore produce an appar- 
ent "floor" in group entropies (C. Norman 2000, 
private comm.; Voit & Bryan 2001). 

Figure 7 plots the entropy at 10% of the Virial 
radius against group X-ray temperature, for the 
simulated groups and the observational analysis of 
PCN99. The group entropy profiles (Figure 5) are 
reasonably well described by power laws. We com- 
pute entropies for Figure 7 by fitting each group's 
profile with a power law, then taking the value 
of the fit at i? = 0.1i?vir- The solid line shows 
the self-similar entropy relation derived using adi- 
abatic simulations (see PCN99), 5'(0.1i?vir) = 
45(r/keV)(/gas/0.06)-2/3/i-4/3]^eV cm'^, with 
/gas = ^b/^m = 0.118 in our case. 

It is clear that the simulated group entropies do 
not follow the self-similar relation, and are instead 
in reasonable agreement with the data (given their 
large error bars). This agreement is somewhat sur- 
prising given that our surface brightness profiles 
do not show the strong mass dependence inferred 
by PCN99, and it once again suggests that the 
canonical beta model may not be an appropriate 
description of poor groups. Our results are also 
in agreement with the analytic model of Bryan 
(2000), which is based on a declining hot frac- 
tion with group mass similar to that predicted by 
our simulation. Thus, it appears that the physi- 
cal mechanism of cooling out low-cntropy gas and 
leaving an effective "entropy floor" (Voit et al. 
2002) is operating in our simulations. Note, how- 
ever, that the simulation does not indicate a hard 
floor at ^ 100 keV cm^, but rathc^r a decline of en- 
tropy with group temperature that is slower than 
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predicted by self-similar models. 



5. Modeling Observed Systems 




-0.5 

log \ (keV) 



In the real universe, it is not possible to directly 
observe the luminosity, temperature, and dark 
matter velocity dispersion of a zero-metallicity in- 
tragroup medium out to the Virial radius. Hence, 
to compare theoretical predictions to observations, 
we must model observational effects in our simu- 
lations. In this section we consider some "real- 
world" effects and assess how they affect the in- 
ferred scaling relations and derived interpreta- 
tions. Our analysis here is simple and intended 
only to indicate the sign and rough magnitude of 
such effects; we defer a detailed side-by-side com- 
parison of simulated and observed groups to future 
work. We will see that there are non-trivial differ- 
ences between the scaling relations of "observed" 
groups and the idealized ones discussed in §3, but 
that none of our fundamental conclusions are al- 
tered as a result. 

5.1. Metallicity 

Since metal lines provide a dominant contribu- 
tion to the total luminosity at temperatures below 
a few keV, they could in principle significantly af- 
fect the predicted luminosity scaling relations. We 
include metals via a fit to the metallicity vs. tem- 
perature data in Figure 1 of Davis, Mulchaey, & 
Mushotzky (1999), namely 



logZ= 1.04 logTx -0.73, 



(6) 



Fig. 7. — Entropy at 0.1i?vir vs. X-ray temper- 
ature of simulated groups. The solid line shows 
the self-similar prediction, and the data points are 
taken from PCN99. 



where Z is the metallicity in solar units, and Tx 
is in kcV. We cap the metallicity at the cluster 
value of 0.3 solar; systems above 1.6 keV are thus 
set to this metallicity. Note that Davis, Mulchaey, 
& Mushotzky (1999) derive a much steeper de- 
pendence of abundance on temperature for sys- 
tems with Tx < 1.5 keV {Z oc T^-^). How- 
ever this relation would produce unrealistically 
low metallicities for our smallest groups, so in- 
stead we use our fit to their entire data set up to 
Tx ~ 3 keV. While there has been some contro- 
versy regarding metallicities derived from ASCA 
data (see, e.g., Buote 2000), the low metallici- 
ties in the Davis, Mulchaey, & Mushotzky (1999) 
sample have now been confirmed in several cases 
with improved data from XMM (Mulchaey, pri- 
vate comm.). Thus, the drop in metallicity below 
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Fig. 8. — Scaling relations as in Figure 1, but 
the X-ray luminosities and temperatures have 
been computed using a metallicity given by equa- 
tion (6). The best- fit relations from Figure 1 have 
been reproduced here as the dotted lines. 



2 keV appears to be real, though it is not straight- 
forward to explain (Mulchaey 2000). 

Figure 8 shows our predicted scaling relations 
after incorporating equation (6) into our calcula- 
tions with the Raymond-Smith code. The rela- 
tions continue to show the same qualitative be- 
havior, and the break values are not changed sig- 
nificantly. The luminosity relations are shallower 
above the break, since the metallicity is mostly 
constant in these larger systems, and metal lines 
provide a greater contribution in smaller (cooler) 
groups. Below the break, the relations are steep- 
ened because the metallicity is dropping with 
group mass. At 1 keV, the luminosity is higher 
by a factor ~ 3 compared to the zero-metallicity 
case (compare the dotted and solid lines in the lu- 
minosity relations), indicating the dominant con- 
tribution of line emission. The Tx — ctdm relation 
is also slightly affected, as higher mass systems 
have lower luminosity-weighted temperatures; this 
is because the relative luminosity weighting of low- 
temperature gas in the system increases with the 
presence of metals in this temperature regime. 
However, these small differences do not qualita- 
tively affect our conclusions drawn from the zero- 
metallicity case (Figure 1). 

While it appears that the metallicity in groups 
is lower than in clusters, this difference is not well 
established. Wc therefore test another model us- 
ing a constant metallicity of 0.3Zq in all groups, 
and we list the resulting scalings in Table 1. This 
calculation shows that if the metallicity is similar 
in groups and clusters, then the Lx — cdm and 
Lx — Tx slopes are shallower in all regimes. Even 
in this case, however, the breaks in the luminosity 
scalings remain quite evident. 

5.2. Surface Brightness Effects 

While current cluster observations can reach 
the Virial radius and sometimes beyond, this is 
rare in groups (Mulchaey 2000). Thus, group ob- 
servations only probe the innermost regions, and 
it is not straightforward to compare their prop- 
erties with clusters where physical quantities are 
measured over the entire system. We investigated 
how this difference affects derived baryon fractions 
in §3.2, and we now consider its effect on scaling 
relations. 

Figure 9 shows the result of applying the surface 
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Fig. 9. Scaling relations as in Figure 1, but the 
X-ray luminosities and temperatures have been 
computed only out to an "observable" fraction of 
the Virial radius as given by equation (1). These 
calculations assume zero metallicity. The best-fit 
relations from Figure 1 have been reproduced here 
as the dotted lines. 



brightness cut in equation (1) to our groups. Note 
that we have reset the metallicity to zero to isolate 
the effects. Comparing to Figure 1 (whose best-fit 
scalings are reproduced as the dotted lines), the lu- 
minosity relations are somewhat steeper, as is ex- 
pected since the surface brightness cut is stronger 
in smaller systems. But the effect is fairly mild, 
even in small groups where the fraction of the 
Virial radius probed is 20% or less. The Tx — <tdm 
relation is essentially unaffected. 

The reason that the surface brightness cut has 
only a mild effect is that most of the group lu- 
minosity comes from the central region (cf. Fig- 
ure 5), so there is minimal impact until the cut- 
off radius becomes a fairly small fraction of the 
Virial radius. Still, one must be cautious about 
interpreting observations where less than 20% of 
the group's Virial radius is being probed, espe- 
cially because the scatter in X-ray siuface bright- 
ness between groups will select out the brightest 
groups of a given mass. Note that we did not 
include scatter in the X-ray extent relation used 
in equation (1), whereas observations show a size- 
able scatter (Mulchaey 2000). Overall, however, 
current observations are sensitive enough to probe 
scaling relations down to fairly small groups rea- 
sonably reliably, assuming that issues of point- 
source contamination and Galactic foreground can 
be overcome. 

5.3. Velocity Dispersion Bias 

In computing Lx — cr and Tx — cr relations in §3, 
we have used the dark matter velocity dispersion 
ctdm- However, the observed velocity dispersion 
is typically computed from galaxies, which may 
be biased relative to the dark matter, and this 
can influence the interpretation of the observed 
correlations. In this section we examine this bias 
and its effect on the observed scaling relations. 

The top panel of Figure 10 shows the group 
mass vs. ctdm • There is a very tight correlation, as 
would be expected for relaxed systems. As before, 
the plot symbols indicate the number of galaxies in 
each group. The bottom panel shows the ratio of 
the galaxy velocity dispersion o-gai to ctdm, where 



(7) 



and (vgai) is the (unweighted) mean velocity of the 
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Fig. 10. — Top panel: Group mass vs. dark matter 
velocity dispersion. Bottom panel: Ratio of veloc- 
ity dispersion calculated from galaxies to that of 
the dark matter. The dotted line shows a linear 
fit to the data points. 



galaxies. There is a weak but discernible trend for 
smaller systems to have their velocity dispersion 
underestimated. The dotted line shows the best 
fit linear relation, which suggests that on average, 
at aDM = 100 km/s, the galaxy velocity disper- 
sion underestimates the true velocity dispersion 
by 15%. Conversely, the galaxy velocity disper- 
sion of large groups is slightly higher than ctdm 
on average. Moreover, there is an increasing scat- 
ter to smaller systems, since there arc only a few 
galaxy tracers with which to estimate the disper- 
sion. This scatter could be larger still in observa- 
tional studies because of projection effects (Tov- 
massian. Yam, & Tiersch 2002). Hence, using 
galaxy velocity dispersions, particularly in systems 
with few identified members, may have an impact 
on the velocity dispersion scaling relations. 

Figure 11 shows the simulated group scaling 
relations Lx — Cgai and Tx — <^ga.i- Lx — Tx is 
not shown since it is (obviously) unchanged from 
Figure 1. The velocity dispersion scalings, on the 
other hand, are substantially altered. Lx — o-gai 
shows a somewhat shallower slope consistent with 
self-similarity, and the fitted slope below the pre- 
vious break value of 180 km/s is now identical to 
the slope above it. Meanwhile, Tx — cTgai shows a 
relation much shallower than obtained using ctdm, 
with Tx oc a'^-Ji^. 

To examine this result in greater detail, we se- 
lected the eleven simulated groups that have ten 
or more member galaxies, then took a subset of 
the A'gai most massive galaxies in each gro\.ip and 
recomputed (Jgsd- lu Figure 12, points with solid 
error bars show the mean and the l-cr scatter of 
the ratio (Jgai/fJDM as a function of the number of 
galaxies iVgai used in the dgai calculation. At high 
A/gai, the ratio is close to unity and the scatter is 
small, but for low TVgai the scatter is large and the 
mean is biased low. If we exclude the most mas- 
sive galaxy from the velocity dispersion calculation 
(points with dotted error bars) , then the mean ra- 
tio remains close to unity even at low iVgai, which 
demonstrates that the bias in cTgai arises mainly 
from the tendency of the most massive, central 
galaxy to move slowly with respect to the group 
ccnter-of-mass velocity. However, if we recalculate 
the Tx — CTgai relation using velocity dispersions 
that exclude the most massive galaxy, we still get 
a fitted slope of « 1.1. This demonstrates that it 
is the random scatter in the galaxy velocity dis- 
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Fig. 11. — Velocity dispersion scaling relations 
with the velocity dispersion computed from the 
galaxies rather than the dark matter. The best-fit 
relations from Figure 1 are shown as the dotted 
lines. 



persion estimates that affects the Tx — cTgai slope, 
not the small systematic bias of cxgai with respect 
to (TDM- Our fits to Lx — CTgai and Tx — o-gai are 
performed in log space without accounting for ran- 
dom errors in the estimation of the true a from 
CTgai. Since downward fluctuations reduce logo- 
more than upward fluctuations increase it, the ef- 
fect of these random errors is to flatten the Lx — cr 
and Tx — cr relations, especially in the low mass 
regime where the errors are larger. 

Noisy velocity dispersion estimates may explain 
some conflicting observational results regarding 
the Tx - cr relation. HPOO flnd Tx oc a°fi, and 
they interpret the higher-than-expected temper- 
atures at low CTgai as evidence for pre-heating. 
However, the velocity dispersions of their smallest 
groups were typically computed from 3-5 mem- 
ber galaxies (compiled from the literature). It is 
not strictly fair to compare their slope with our 
Tx — CTgai slope, since our galaxy samples have dif- 
ferent "completeness limits", but it is illustrative 
that the slopes are similar. Conversely, the MZ98 
sample includes only groups with over 20 spectro- 
scopic members, and they flnd Tx oc cf'^^i^'^'^, in 
agreement with cluster samples. 

In the currently relevant observational regime 
(ct > 200 km/s), the Lx — cr relation is not signif- 
icantly affected by the scatter in cTgai (Figure 11, 
top). HPOO find Lx cx CTggl, consistent with self- 
similarity, but they attribute this to coinciden- 
tally canceling reductions in Lx owing to pre- 
heating, and in CTgai owing to a noisy estimation 
from a small number of galaxies. Our analysis 
shows that, at least for groups of this size, the 
Lx — CTgai relation is unaffected within current ob- 
servational errors. In support of this idea, MZ98 
and Zimer, Zabludoff & Mulchaey (2001) also 
find Lx cc CT^.^^^'* using many more galaxies per 
group. Tovmassian, Yam, & Ticrsch (2002) argue 
that the velocity dispersion bias arises because of 
geometrical projection effects in 1-D dispersion es- 
timates, and they also find Lx oc ct^^^j when prop- 
erly accounting for such effects. 

In summary, velocity dispersion bias can have 
a significant effect on the Tx — ct relation, but 
only a mild effect on Lx — cr, at least above the 
break. Observations that have sufficient numbers 
of galaxies per group (> 10) indicate little evi- 
dence for the anomalously hot intragroup gas in 
poor groups that would be expected in a pre- 
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heating scenario. 



5.4. Comparison to Observations 
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Fig. 12. — Velocity dispersion calculated from the 
A/gai most massive galaxies within a group, di- 
vided by the velocity dispersion of the dark mat- 
ter. Points with solid error bars show the calcula- 
tion using equation (7), while points with dotted 
error bars indicate the effect of excluding the most 
massive galaxy before computing o-gai- Results are 
computed only for the 11 simulated groups that 
have ten or more member galaxies. The error bars 
represent the 1 — a scatter of CTgai values within 
each bin; the error on the mean of tTgai/uDM is 
smaller (e.g., for points at N^ai < 10, it is lower 
by a factor of v^). 



Though there are still many fundamental dif- 
ferences between the way our simulations and ob- 
servations are analyzed, as well as some inherent 
modeling difficulties in our simulation, it is never- 
theless instructive to compare our scaling relation 
predictions with available observations. To model 
observations as closely as we can given current 
constraints, we apply the metallicity relation as in 
equation (6), we include surface brightness effects 
as in equation (1), and we use the dark matter 
velocity dispersion. The last choice is motivated 
by the fact that the latest group samples have 
deep imaging that permits identification of 10-50 
group members (e.g. Zimer, Zabludoff & Mulchaey 
2001), which should be sufficient to trace the true 
mass of the group much like ctdm, though still with 
somewhat larger scatter. 

Figure 13 shows the resulting scaling relations, 
and the fits are also listed in Table 1. The Lx — 
f DM relation slopes above and below the break are 
close to those in a compilation of clusters, groups 
and galaxies by Mahdavi fc Geller (2001), who 
found Lx oc cr^'^ down to a = 350 km/s and 
Lx oc below that down to galaxy scales. Our 
break occurs at a lower a than theirs, but the ob- 
served relation is poorly constrained in the break 
region, and the data sample is inhomogeneous, so 
the discrepancy may not be too serious. 

The slope of the Lx — Tx relation above the 
break is in good agreement with MZ98 and Zimer, 
Zabludoff & Mulchaey (2001), and it is also 
consistent with the observed relation in clusters 
(White, Jones, & Forman 1997). However, it 
is in poor agreement with samples by HPOO and 
Xue & Wu (2000), who find Lx oc T^^ and 
Lx oc T^^. Those results are in better agree- 
ment with the slope helow our break at 0.8 keV. 
These samples have Tx ~ 0.4—1.6 keV, thus their 
steep slopes suggest a higher break temperature, 
perhaps ~ 1 — 1.5 keV. As with our Lx — o rela- 
tion, our break appears to occur at a mass that is 
slightly too low. 

Our Tx — o relation shows some deviation from 
self-similarity (Ty oc cr^'^^), although as we saw 

in §5.1, this deviation arises mostly from metal 
line emission weighted into the temperature de- 
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Fig. 13. Scaling relations as in Figure 1, but the 
X-ray luminosities and temperatures have been 
computed using a metallicity given by equation (6) 
and a surface brightness cut given by equation (1). 
This represents our attempt to most closely match 
current observations, which are discussed in the 
text. 



termination, not from any physical process that 
actually raises the gas temperature. Metal line 
cooling is another way in which the Tx — <J rela- 
tion can mimic excess temperature in Tx ^ 1 keV 
groups, in addition to the effects of scatter in a 
estimates from small numbers of galaxies (§5.3). 
Observations suggest that the Tx — cr relation is 
fairly close to self-similar from clusters down to 
groups (Mulchaey 2000), in agreement with our 
predictions. 

So far we have not discussed the amplitudes 
of the various scaling relations, only the slopes. 
Our Lx — cr relations amplitude is in very good 
agreement with observations. MZ98, HPOO, and 
Zimer, Zabludoff & Mulchaey (2001) all find 
Lx = 10''2-5-10^2-''erg s'^ at a = 300 kni/s where 
the observations are relatively reliable, while we 
obtain Lx = 10^^'^erg s~^. This is encouraging, 
given that this (large-volume) simulation cools an 
excessive fraction of baryons (see Figure 2 and 
the accompanying discussion), and we have made 
large corrections to the luminosities based upon 
our two-phase decoupling (see §2). If our predic- 
tions were highly sensitive to these effects, it would 
be quite a coincidence that they conspired to bring 
our luminosity-mass relation into such good agree- 
ment with observations. 

Conversely, the amplitudes of the temperature 
relations are in poor agreement with observations. 
As discussed in §4.3, our luminosity-weighted tem- 
peratures are higher than observed by a factor of 
« 1.5 — 2 at a given cr. The lower factor comes 
from considering (3spec as observed in Chandra 
data of large clusters by Allen, Schmidt, & Fabian 
(2001) (with the caveat that our systems are much 
smaller than any observed by those authors); the 
higher factor comes from comparing the amplitude 
of our Tx — <y relation with iiOS'^T observations. 
For example, we predic:t a temperature of 1.7 keV 
at CT = 300 km/s, while MZ98 and HPOO find 0.8 
and 0.9 keV, respectively, for such systems. It re- 
mains to be seen whether Chandra observations of 
these groups will yield higher temperatures as has 
happened for some clusters. 

In the Lx — Tx relation, this overestimation 
of temperatures translates to our predicted lu- 
minosities being too low by a factor ^3 — 8 
at fixed Tx (for Lx oc T|-). For example, at 
Tx = 1 keV, MZ98 find Lx = W'^'^ '^erg s'^, HPOO 
find Lx = lO^^ergs"^, and the extrapolation of 
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the White, Jones, & Forman (1997) relation yields 
Lx = 10**^ ''erg s~^; our simulated groups yield 
Lx = 10^^'^erg s~^ at this temperature. As dis- 
cussed in §4.3, the ad hoc correction of lowering the 
temperature of gas within 0.2i?vir in accordance 
with temperature profiles seen by HPOO results in 
lowering our groups' Tx by ^70%, and such a re- 
duction (regardless of what physical phenomenon 
is responsible for it) brings both our Tx — o and 
Lx — Tx relations into the observed range. 

In summary, our simulation seems to do a good 
job of reproducing the slopes of observed scal- 
ing relations, but the breaks in the Lx — cr and 
Lx — Tx relations appear to occur at a group 
mass that is somewhat too small. This discrep- 
ancy might be resolved by improvements in the 
numerics or astrophysical modeling, or it could in- 
dicate that some degree of pre-heating is still re- 
quired, though much less than in a model with adi- 
abatic gas dynamics. The amplitude of our Lx—cr 
relation agrees reasonably well with observations, 
suggesting that our groups' luminosities are be- 
ing accurately modeled. However, our simulated 
X-ray temperatures appear to be discrepant with 
observations, being too high by a factor of 1.5 — 2 
at a given mass. This discrepancy manifests itself 
in both the Lx — Tx relation and the Tx — (J rela- 
tion. The lower temperatures of observed groups 
appear to be connected to cooler central regions 
that are not reproduced by our simulation. We 
suspect that the cause of this discrepancy is astro- 
physical, perhaps our use of zero-metallicity cool- 
ing curves during dynamical evolution, or perhaps 
the absence of a significant physical process such 
as AGN heating or conduction. 

6. Conclusions 

Using a large-scale hydrodynamic cosmological 
simulation that is able to resolve sub-L* galaxies, 
we have investigated the X-ray scaling relations 
for galaxy groups with velocity dispersions ranging 
from 100 — 550 km/s. In particular, our goal was 
to assess whether these scaling relations are con- 
sistent with those predicted by simple self-similar 
models and adiabatic simulations. Our main find- 
ings are: 

1. Our luminosity scaling relations Lx — Tx 
and Lx — cr do not follow the self-similar 
form at any point in our mass range. The 



simulated relations are always steeper than 
the self-similar predictions, with lower lumi- 
nosities at smaller masses. 

2. These scaling relations show a "break" 
around 180 km/s (0.7 keV), below which 
the relations steepen even further. 

3. The Tx — cr relation shows no break, and its 
slope is reasonably consistent with the self- 
similar model, indicating that these groups 
are Virialized and in hydrostatic equilib- 
rium. 

4. In our simulations, self-similarity is broken 
mainly as a result of the hot gas fraction 
falling with group size, from « 50% of the 
baryonic mass at a = 500 km/s to « 20% 
at (7 = 100 km/s. There is a secondary ef- 
fect from the fact that the hot gas clumping 
factor (mostly in the groups' central regions) 
drops by a factor of roughly three over this 
mass range. 

5. Our large-volume simulation cools too many 
baryons, resulting in a typical hot frac- 
tion that is lower than observed by 40%, 
once X-ray surface brightness thresholds are 
taken into account. Our high-resolution sim- 
ulation predicts more hot gas at a given 
group mass, and is in good agreement with 
observations. Thus, the amplitude of the 
/hot — cr relation is somewhat sensitive to 
resolution, though not in the sense that one 
might naively expect. The trend of falling 
hot gas fraction towards lower group mass, 
which arises in both simulations and ac- 
counts for most of the departure from self- 
similar scaling relations, is consistent with 
available observations. 

6. Our two simulations, which differ in mass 
resolution by a factor of eight, predict very 
similar X-ray luminosities at fixed group 
mass. To the extent that we can test them 
with these simulations, the departures of the 
Lx — cr and Tx — cr relations from self-similar 
scaling are not artifacts of numerical resolu- 
tion: the reduced luminosities of low mass 
groups are a result of physical processes that 
are modeled consistently between the two 
simulations. 
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7. Surface brightness profiles offer an excellent 
diagnostic for constraining gas physics in 
groups, in particular for distinguishing cool- 
ing from pre-heating as the primary cause of 
departure from self-similar scalings, but cur- 
rent observations are too uncertain to give 
clear guidance. Our groups' median [3^ = 
0.64 is intermediate between various obser- 
vations, and we predict no significant trend 
of Pat with group size. Our simulations do 
reproduce the observed entropy-temperature 
relation (PCN99), though we do not predict 
a hard floor but merely a shallower decline 
of entropy with group mass than predicted 
in self-similar models. 

8. Accounting for the effects of metallicity or 
for observational biases due to X-ray surface 
brightness thresholds does not qualitatively 
alter the conclusions listed above. 

9. Noisy estimates of group velocity disper- 
sions have a significant effect on the derived 
Tx ~ cr relation. The most massive galaxy 
in a group tends to move close to the group 
center-of-mass velocity, so its inclusion in 
a velocity dispersion estimate biases a low 
when the number of galaxies is small. Even 
more important is the scatter in ct when the 
number of tracers is small, which biases a 
simple fit of Tx — o. Our simulations sug- 
gest that > 10 galaxies are needed to obtain 
velocity dispersions that are precise enough 
for Tx — cr determination. Noisy a estimates 
may account for some of the discrepancies 
between apparently conflicting observational 
analyses. The Lx — cr relation can also be 
biased by noisy a estimates, but this effect 
does not seem to be important in the mass 
range probed by current data. 

10. Accounting for observational biases, our 
scaling relation slopes are roughly in agree- 
ment with observations, though our break 
appears at somewhat too low a group mass. 
With our best attempts to mimic current 
observational approaches, our simulation 
predicts the following scaling relations for 
galaxy groups with a < 500 km/s: 

• Lx oc a^'* for a > 180 km/s, and 
Lx cx ct" '^ for a < 180 km/s. 



• Lx oc r|-6 for Tx > 0.7 keV, and 
Lx oc r|-8 for Tx < 0.7 keV. 

• Tx oc a^'"^^ when > 10 galaxies used to 
estimate a, and 

Tx oc cr~^ when 3-5 galaxies are used 
for the smallest groups. 

The relations with various assumptions 
about metallicity, surface brightness thresh- 
olds, and velocity dispersion estimates are 
summarized in Table 1. 

11. While the amplitude of our Lx — cr relation 
amplitude is in good agreement with obser- 
vations, the amplitudes of our Tx — cr and 
Lx — Tx relations show significant discrep- 
ancies, and our median value of /Jgpoc = 0.47 
is lower than observational estimates. We 
conclude that the luminosity-weighted tem- 
peratures of our simulated groups are too 
high by a factor of 1.5 — 2. This discrepancy 
appears to originate in the central regions, 
with r < 0.2_Rvir, where observed group 
temperature profiles drop towards the cen- 
ter while our simulated profiles are flat or 
rising. An ad hoc introduction of a cool cen- 
tral region brings the amplitudes of both the 
Tx—cr and Lx—Tx relations into reasonable 
agreement with observations. We suspect 
that the ultimate origin of these discrepan- 
cies is missing physics in our simulation, such 
as AGN heating of the intracluster gas, but 
we cannot rule out numerical effects. 

These results can be taken on several levels. 
At their most basic and most robust level, they 
indicate that groups cannot be treated as self- 
similar ly scaled-down versions of clusters, since 
radiative cooling plays an increasingly important 
role in smaller systems. Although the cooling time 
in the outskirts of groups and clusters is longer 
than the Hubble time today, cooling times were 
shorter in the subsystems that merged to make 
up the final group, and even today cooling can 
be significant out to a non-trivial fraction of the 
Virial radius. Thus, the conventional approach of 
modeling X-ray properties with purely adiabatic 
physics must be applied with caution, especially 
as one moves from the mass regime of rich clus- 
ters to that of poor clusters and groups. Since 
cooling is a process that is known to occur, more 
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exotic processes like pre-heating or entropy injec- 
tion should be examined in the context of models 
that already incorporate cooling. Unfortunately, 
this requirement increases the complexity and the 
uncertainty of the calculations, whether they are 
analytic or numerical. 

At a second level, our results show that this sim- 
ulation with cooling qualitatively reproduces the 
trends seen in observed scaling relations, namely 
a reduced luminosity of lower mass systems (rel- 
ative to self-similar predictions), and a break be- 
low which the luminosity relations steepen further. 
While we have used numerical methods to obtain 
this result, it may be explained by a simple phys- 
ical picture based on hierarchical growth of struc- 
ture. At any epoch, the progenitor of a group 
will have a lower Virial temperature than the pro- 
genitor of a cluster. Hence, integrated over the 
formation history of these objects, a larger frac- 
tion of baryons will cool in the group environment 
than in the cluster environment. This results in 
the trends seen in the simulations and in the ob- 
servations. Put another way, it is well known that 
in clusters most baryons are hot, while in galaxies 
most baryons are cold; our results simply suggest 
that the transition between these regimes occurs 
in galaxy groups around 1 keV. 

At a third, more uncertain level, our results 
may be taken as evidence that a ACDM cosmo- 
logical model with standard gas dynamics and ra- 
diative cooling can quantitatively explain observed 
group scaling relations, without strong pre-heating 
or substantial entropy injection. The uncertainty 
arises because our simulation does not reproduce 
the observations in their entirety: the forms of the 
predicted scalings are about right, but the mass 
scale of the breaks is somewhat too low, the cold 
gas fractions are somewhat resolution-dependent, 
and at a given a the luminosity-weighted temper- 
atures are too high. Presently, we do not know the 
extent to which these discrepancies reflect numer- 
ical inaccuracies, missing astrophysical processes, 
incorrect cosmological parameters, or errors in the 
observational inferences themselves. We therefore 
do not know whether solving these problems will 
reduce the role of cooling to the point that it no 
longer explains the observed form of the scaling re- 
lations. Definitive measurements of resolved sur- 
face brightness and gas density profiles would be 
a powerful diagnostic for the relative importance 



of radiative cooling and non-gravitational heating 
in accounting for observed scaling relations. 

Our simulation demonstrates that cooling can 
reproduce many of the qualitative features in scal- 
ing relations that are often quoted as evidence for 
pre-heating. Our conclusion on this point agrees 
with those drawn by Bryan (2000) and Voit & 
Bryan (2001) on the basis of analytic models and 
by Muanwong et al. (2001) on the basis of nu- 
merical simulations. New in this paper are the 
use of relatively high-resolution simulations well 
suited to the group mass regime and a more de- 
tailed consideration of observational issues than 
these earlier studies. In agreement with Bower et 
al. (2001), we find that cooling seems to "kick 
in" at a mass scale slightly smaller than required 
by observations. This discrepancy of mass scales 
may indicate that some non-gravitational heating 
is still required. Certainly high-redshift galaxies 
are observed to produce strong winds (Pettini et 
al. 2001), and AGN activity is sometimes seen in 
the centers of groups, so potential sources of non- 
gravitational heating do exist. Some recent mod- 
els of cluster "cooling flows" incorporate an addi- 
tional heat source that serves to maintain the inner 
gas at an intermediate temperature, suppressing 
the spectral signatures of ~ 1 keV gas that would 
otherwise be expected (David et al. 2001; Fabian 
et al. 2001). Hence, it is plausible that baryons 
in groups have experienced some non-gravitational 
heating, but the amount needed to reconcile mod- 
els with observations may be much less than previ- 
ously proposed, thus alleviating the "ICM energy 
crisis." 

Fortunately, this field appears poised for break- 
throughs on both observational and theoretical 
fronts. XMM and Chandra will allow studies 
of groups with greater sensitivity, spatial resolu- 
tion, and spectral resolution than previously pos- 
sible. Constraints on optical properties of groups 
continue to improve with deeper surveys. On 
the modeling side, much effort has gone into ex- 
ploring solutions for the "decoupling problem" 
in two-phase regimes, with new algorithms be- 
ing proposed that may alleviate the problem 
in a self-consistent fashion (e.g. Thacker et al. 
2000; Springel & Hernquist 2001). Furthermore, 
Moore's Law improvement in computing technol- 
ogy continues unabated, enabling ever larger and 
higher resolution simulations. Advances on both 
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fronts will soon enable reliable sidc-by-side com- 
parisons of simulations and observations, where 
the simulations are analyzed using exactly the 
same techniques as the data. This approach 
promises to yield an in-depth understanding of 
the physical processes governing the formation 
of galaxies in their most ubiquitous environment, 
groups. 
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